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In this paper, we consider a time-headway model, introduced in Ref. |4§, for buses on a bus route. 
By including a simple no-passing rule, we are able to enumerate and study the unstable modes of a 
homogeneous system. We then discuss the application of the model to realistic scenarios, showing 
that the range of reasonable parameter values is more restricted than one might imagine. We end 
by showing that strict stability in a homogeneous bus route requires careful monitoring by each bus 
of the bus in front of it, but in many cases this is unnecessary because the time it takes for the 
instability to appear is longer than a bus would normally spend on a route. 



I. INTRODUCTION 

While there has been much interest in the study of automobile traffic pj , there have been few corresponding studies 
of buses 

H H II The dynamics of a bus route, while having some similarities with that of automobile traffic, 
differs due to the added interaction of buses with passengers at designated bus stops. A good reason for studying the 
dynamics of bus routes is that they are so often unstable. Buses are initially spaced at regular intervals. However, 
if one bus is delayed for some reason, it will find a larger number of passengers waiting for it at subsequent stops, 
delaying it further. Meanwhile, the bus following finds fewer passengers waiting for it, allowing it to go faster until 
eventually it meets up with the delayed bus. Clusters of three, four, or more buses have been known to form in this 
manner, resulting in slower service. 

In references [2j and |J|, Nagatani presents a time-headway model for buses. Using linear stability analysis, he is 
able to determine the range of parameters over which the homogeneous solution (i.e., with buses spaced evenly apart) 
is unstable. In this paper, we make a more thorough investigation of Nagatani's model. We demonstrate the existence 
of three types of phase diagrams, in which the behavior of the bus system is divided into four separate categories. We 
conclude with a discussion of how this model may be applied to real- world situations, and the limitations imposed by 
practical considerations. 

II. MODEL 

We consider the following model, introduced in of buses on a bus route (Fig. Bus stops are labelled by 
s = 1, 2, . . . where stops s and s + 1 are a distance L apart. There are J buses, j = I, . . . , J, which travel from stop 
to stop, with bus j = 1 in the lead and bus j = J in the rear. Every bus visits every stop, and buses do not pass one 
another. The time at which bus j arrives at stop s is tj jS , which is given by the recursive relation 

tj, s - *i,s-i = <Vy<% s _i + — , (1) 

V J;S-l 

where 

$tj.s — tj s tj — l.s (2) 

is the time-headway, the time gap in front of bus j at stop s. The first term on the right-hand side of Eq.^is the time 
it takes for passengers to board the bus at stop s — 1. The parameter A is the rate at which passengers arrive at a bus 
stop; XStj s _i is the number of passengers that have arrived at stop s — I since the previous bus left. The parameter 
7 is the time it takes each passenger to board the bus, so \j5tj tS -i is the amount of time needed to board all of the 
passengers. For convenience, we introduce the dimensionless parameter /i = A7, which we call the passenger rate. For 
simplicity we ignore the passengers getting off of the bus; we will assume that it takes less time for the passengers to 
get off than it does to get on and pay the fare. 



2 



s-3 



s-2 



s-1 



s + 1 



s + 2 



s + 3 



L 



L 



L 



bus 
j+1 

[□on 



bus 
j 



bus 
jl 



1 



Atj+l 



At 



FIG. 1: Schematic illustration of the model. 



The second term in Eq.^is the time it takes for bus j to travel from stop s — 1 to stop s, where Vj^—i is the average 
velocity of the bus between stops. If this velocity is constant, then the tendency for buses to bunch together, as 
described in the introduction, has no counterweight, and a steady flow of buses will always be unstable (unless there 
are no passengers). It is reasonable to assume, however, that a bus driver will try to prevent bunching by slowing 
down when the gap between his bus and the next is too small. One can model this by writing the average speed Vj, s 
as a function V(Stj S ) of the gap between his bus and the bus in front of him: 



V(5t) = v min + (v max - v min ) 



tanhw(5i — t c ) + tanhwt c 
1 + tanh u)t c 



(3) 



The hyperbolic tangent factor acts as a spread-out step function, centered at t c with a width proportional to 1/oj. 
The parameter t> max is the speed a free bus (i.e. one that is alone on the route) would travel. On the other hand, i> m i n 
is the speed a bus travels if it has completely caught up with the bus in front of it. For example, if v m i n = then a 
bus which has caught up with the bus in front of it will stop and wait until its predecessor has cleared the next stop 
before proceeding. 

In what follows, it is convenient to work with the time headways Stj jS , rather than the arrival times tj jS . It is also 
convenient to rewrite our expressions in terms of dimensionless quantities. In doing so, we find that there are four 
significant parameters, not including initial conditions. The first such parameter is the passenger rate fi. The other 
three are a = Lw/f ma x, P = v m ya/v mlCK , and e = 1 — tanhwt c (which will typically be small). We will also consider the 
dimensionless variable Atj tS = uj5tj jS and the dimensionless velocity function V(At) = V(At)/v max . Our evolution 
equation Eq. ^ now reads 



where 



At,- s - A< 7 -, s _i = a 



V(At) =(3 + 



1 



1 



- /i[Ai, )S _i - Aij_i )S _i], 



V(Atj, a -i) V(A^-m-i). 
(1 - /3)e tanh At /3(1 - tanh At) + e tanh At 



1 - (1 - e) tanh At (1 - tanh At) + e tanh At 



(4) 



(5) 



III. STABILITY ANALYSIS 



We are interested in the stability of a homogeneous flow of buses, with Atj jS = Atj$ = Ato. One can easily 
verify that this is a solution to Eq. 0] Starting with a small perturbation to the initial homogeneous solution: 
Atj s = Ato + Uj.s, where yj a is small. To first order, Eq. 0] becomes 

Vj,s - Vj,s-i = [Vj+i.s-1 - yj,s-i][F(At ) - fi\, (6) 
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FIG. 2: Phase diagrams indicating the regions of stability of Eq. |S| for three representative values of the parameters. The 
shaded area is the region which satisfies Eq. [HI 



where we have introduced the convenient abbreviation 

V'{At Q ) a(l - f3)e(l - tanh 2 Ai ) 



F(Aio) = a- 



' V(At ) 2 [/3(1 - tanh At ) + e tanh At ] 2 ' 
It can be shown Q that the perturbation is stable if 

F(Ato) - 1< m < F(At ). 



(7) 



(8) 



From Eq. [H] we can construct a phase diagram (Fig. [2J for the stability of an initially homogeneous bus route, 
based on the initial spacing Atn and the passenger rate /x. The stable region in phase space is bounded by the curves 
/i = F(Ato) and /i = F(At ) — 1. Because of the added constraint that /i > 0, there are different phase diagrams 
depending on whether F(A£n) — 1 is ever positive (Fig. [30 or not (Fig. [3d). The curve F(At) has a maximum value 
of 



(1-/3) e 



so the phase diagram resembles Fig. [3, whenever 



a > 



2/3 -e 2/3 



1-/3 1-/3 



(9) 



(10) 



A third phase diagram, Fig. [3, arises when Umin — 0; as it is in figure 3 of Ref. Q (although apparently not in figure 
8 of the same reference, which may account for the discrepancy between those two phase diagrams.) 



IV. SIMULATION 



To study the ways in which the system becomes unstable, we evaluate Eq. Q] iteratively in s. Our initial condition 

is 

At jt0 = At +0.1^, (11) 

where rj are random numbers chosen between —1 and 1. For each combination of initial headway Ai and passenger 
rate /i, we run the simulation until either a) we reach stop s = 5000, or b) one or more of the bus headways exceeds 
At = 1000 (in which case the system has become unphysical). 
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FIG. 3: Phase diagram for bus systems with (a) periodic and (b) fixed boundary conditions, where a = 1, j3 = 1/4, and 
e = 0.036. The horizontal axis is the initial time headway Ato, while the vertical axis is the passenger rate [i. Stable runs 
are marked by circles (O) an d exploding runs (§ HVBfl by exes (x). Oscillatory solutions (§ HVDfl are marked by 
squares (□); diamonds (o) mark runs which started like oscillatory solutions but ended up flat. Slowed solutions f S II V" (Jl with 
clusters are marked by upward-pointing triangles (A) and slowed solutions without clusters by downward-pointing triangles 
(v)- The grey shading shows the region where F(Ato) — 1 < fj. < F(Ato). 



In this paper we consider two different boundary conditions. The first is periodic in the bus number j; so for 
example Atx lS = ti lS — tj tS . This is convenient numerically, and it creates translational symmetry, but it is hard to 
construct a physical model which begins with this characteristic. We also consider a fixed boundary condition, where 
Ati jS = AtQ' Since the velocity of a bus depends entirely on At, this corresponds to a scenario where the initial bus 
(j = 1) moves at a constant speed V(A£o)-@] 

The structure of the model requires that buses not pass one another; however, there is nothing in Eq. 0|to prevent 
the headways At from becoming negative. To fix this, we add to our simulation the rule that any Atj jS < is replaced 
by Atj^ s = 0. This corresponds to a situation where drivers are forbidden (or unable due to road conditions) to pass 
one another. 

Figure |3| shows the results of our simulation runs for a typical set of parameters (a = 1, f3 = 1/4, e = 1 — tanh2 = 
0.036), using both types of boundary conditions. In both cases, the phase space is divided into four regions, corre- 
sponding to four types of runs. 



A. Stable Runs 



Most of the runs within the stable region, as defined by Eq. [SI remain homogeneous. In the periodic case, the 
initial fluctuations in Atj settle into a small precessing sinusoidal perturbation which decays exponentially with time 
(Fig. UJ. In the fixed case, the system quickly locks onto the constant solution At = Ato with no fluctuations. 



B. Explosive Runs 



Most of the runs lying above the stable region quickly develop an unphysical instability. This takes the form seen in 
Fig. [S] independent of boundary condition: those headways lying above the mean increase exponentially, while those 
lying below decrease steadily until they reach zero. An observer stationed at a stop far down the line will see clusters 
of buses arriving after long waits; far enough down the line, these waits become astronomical, which is absurd. Clearly 
this model is insufficient to deal with these runs at long times. 
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FIG. 4: A plot of max(Atj ) — min(Atj ) versus iteration step s shows that the periodic system (/x = 0.8, Ato = 1.5) is converging 
exponentially to the homogeneous solution. The same system with the fixed boundary condition converges much more quickly. 
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FIG. 5: An extreme example of an explosive run, with /i = 1.9 and Ato = 2.5. By stop s — 8 there are buses which are already 
1000 time units apart (where one time unit is the time it takes for a free bus to travel from one stop to the next). 
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C. Slowed Runs 
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FIG. 6: An example of a slowed run, where /i = 0.95 and Ato = 0.2. 

To the left of the stable region are runs which develop an alternative stable solution, as seen in Fig. [5] In the case 
of the fixed boundary condition, these runs have two things in common. The first, indicated by the vanishing of one 
or more headways, is the appearance of clusters: two or more buses which travel along as a single unit. The second is 
that the units, whether single buses or clusters, are homogeneously spaced, but with a spacing that is larger than the 
initial spacing Ato- ft should be pointed out that the solutions shown here are stationary; the clusters and spacings, 
after they form, do not change. 

In analytic terms, these states arc of the form Atj jS = Atj — rrj, where rj is either or 1, and r > Ato is a 
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FIG. 7: The solution to Ea. 1131 for the set of parameters a = 1, /3 = 1/4, e = 0.036. The x marks the solution corresponding 
to Fig. El m simulation, the spacing r between slowed buses always comes from the lower branch of the curve. 
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FIG. 8: An example of the oscillatory condition, with /i = 0.1, Ato = 1.0, and fixed boundary conditions. This plot shows how 
long each of three consecutive buses arrive at a stop s after the initial bus j = 1 arrived. Notice that the middle bus is bunched 
first with the bus preceding it, then the bus following, and so forth. 



constant. It is straightforward to show that this is a solution to Eq. 0] 



= a 



+ V>T[ r i -rj-i] 



(12) 



When rj = r_y_i, this equation is satisfied trivially. Otherwise, the equation takes the form 



_ a ( 1 1 



(13) 



which we can solve numerically for r (Fig. EJ. For a given passenger rate /i, these spacings r correspond precisely 
with those seen in simulation. Furthermore, for high enough passenger rates — \x > 1.199 for this set of parameters — 
Eg. 1131 has no real solutions, which explains the cut-off in Fig. [3] between the slowed and explosive regimes. 

In the case of the periodic boundary condition there are cases where the clusters eventually break up, leaving a 
system of buses which are equally spaced, but with the larger spacing predicted by Eq. 1131 These runs are marked 
by downward-pointing triangles (y) in Fig. [31 



D. Oscillatory Runs 

In the case of a run lying below the stability region in phase space, the first term in Eq. ^ which is meant to 
resist the tendency for buses to cluster, becomes too large. This leads to overreaction, so that two buses which 
arrive too close together at one stop are too far apart at the next. The resulting behavior may be compared to a 
system of underdamped oscillators. Figure |S] shows the resulting behavior. For the periodic boundary condition, 
these oscillations decay as a power law(Fig. |5J), but at so slow a rate as to be practically permanent. With the fixed 
boundary condition, the earliest buses (i.e., those with the lowest j) shed the oscillating behavior after only a few 
stops, resuming a homogeneous configuration; with more iterations, more buses join the homogeneous regime. In 
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FIG. 9: The decay of an oscillatory run with periodic boundary conditions, [i — 0.2, Ato — 1.2. Both axes are logarithmic. 
The power-law decay is too small to have an appreciable effect on the behavior of the buses. 
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FIG. 10: The headways of an oscillatory system with fixed boundary conditions, fi — 0.1, Ato = 1-0, at four different bus stops. 
The bifurcation of each line indicates the presence of oscillatory behavior, which we show in full for the s = 200 case. 
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some cases, such as in Fig.^J the system reaches a steady state with the oscillations still dominating the later buses. 
In other cases, however, the system becomes completely homogeneous, as the bifurcation point seen in Fig. I1UI slips 
off the right side of the graph. This later effect may be due to the finite number of buses. 



V. DISCUSSION 



In our simulations we have considered a large range of values for /x and Aio. However, these parameters should be 
limited by a couple of practical concerns. 

The passenger rate /i is defined as the product of the number of passengers that arrive at a stop per unit time, 
and the time it takes a single person to board the bus. Put another way, it is the ratio of the number of people 
that arrive at a stop to the number of people who can board the bus in the same amount of time. If this number is 
greater than 1, then passengers arrive at a stop faster than the bus can take them on, and the bus should never be 
able to leave the stop. Since < fi < 1, only one of the two inequalities in Eq. [HI is meaningful for a given value of 
the parameters (including Aio), depending on whether _F(Aio) is larger or smaller than 1. This suggests that if one 
wanted to maximize the area of the stability region in phase space, one would do well to make sure that the lower 
stability curve F(At ) — 1 just grazes zero, or that a(l — f3) = 2(3 — e according to Eq. [3J 

Another practical consideration puts a limit on the value of Ai . Typically, buses are spaced far enough apart so 
that the first bus will reach the first stop before the second bus is allowed to leave, particularly if the stops are spaced 
fairly close together. This is described by the inequality 

L (1 - tanh Ai ) + etanh Ai . . 

Ai > ^ = aA- — -- > a. (14 

V(At ) /3(1 -tanh Ai ) + etanh Ai 

For the parameters we have been studying, 

1 -0.964 tanh Ai * 

Ai > r-r- => Ai > 1.82. (15) 

0.25 - 0.214 tanh Ai 

This cuts out much of the interesting part of Fig |3 including the slowed runs and almost all of the underdamped 
solutions. In our discrete model, the basic iteration step is the bus stop; bus drivers are allowed to change their speed 
at the bus stops and nowhere else. If the buses are several stops apart, then they have enough time to react to one 
another. Otherwise, unusual situations such as the slowing case or the underdamped case may arise. 

Finally, we consider the relationship between Stg and i c , which is how close a bus will come to the bus in front of 
it without slowing down. Consider the ratio 

1 — tanh Aio 1 — tanh ujStn 

(16) 



1 — tanh Lut r 



We can rewrite F(Aio) in terms of this ratio: 



F(At ) = ^-/ Vti~-V 
If er = 1 — tanh Ai <C 1 (which it will be if Eq. Elis valid, since 1 — tanh 1.82 = 0.05), then 

F(Ato)«a(l-/3)^-^. (18) 

Now let us consider what values our parameters might take in real life. A typical urban bus route might have 
L = 0.5 km, w max = 50 km/hr = | km/min, v m [ n = 15 km/hr = j km/min, and u = 1 /min; thus a — 0.6 and (3 = 0.3. 
A bus which runs every 10 minutes might take on two passengers at every stop, so A = 0.2 people per minute. If it 
takes 7 = 3s for a person to board a bus, then ^ = A7 = 0.01. 

Consider a scenario where bus drivers only react to what they see; that is, a driver will only slow down if she 
can see the next bus in front of her. It takes a free bus L/v max = 0.6 minutes to travel from one stop to the 
next, so a reasonable value for the amount of warning a bus driver has is on the order of t c = 1 min. Typical bus 
routes have buses which are spaced much farther apart, perhaps every Aio = 10 min or more. In this scenario, 
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FIG. 11: This shows the number of stops, for a given passenger rate n, that the bus route will run before seeing a 1-minutc, 
5-minute, 10-minute, or 15-minute deviation in the initial homogeneous state. The parameters used in this simulation are 
a = 0.6, /3 = 0.3, e = 1 — tanh 1, and Sto = 60 min. A power-law fit to all four lines shows that they all go as f^- - 965 ^ 005 _ 



r = (1 - tanh 10)/(1- tanh 1) = 10~ 8 , so F(At ) « 10~ 8 . Since /i = 10~ 2 , the stability condition Eq.|H|is very clearly 
violated. For F(Ato) to reach a high enough value to create stability, we need in general for the ratio r to be closer 
to 1. F(r) takes its maximum value when r = 1/(2/3), in which case F(r) — 1/(8(3) = 0.4, which is easily larger 
than /j, in this example. Notice that, for values of t c and Sto greater than 1 min, r f» e 2 ( tc ~ 5t °', so f or each minute's 
difference between t c and Sto, r is increased or decreased by a factor of ten. It would seem that to maintain a stable 
homogeneous bus route, drivers must be reacting to the bus in front of them even from the very beginning, and can 
only ignore the leading bus if they have gotten far enough behind (in which case, of course, the proper solution is to 
go as quickly as possible). 

This is quite a stringent requirement for stability, and explains why it is so common to see clusters of buses in large 
cities. It does not seem likely, however, that this would be the case for less frequent bus service, such as when buses 
run once per hour. A driver on such a route does not typically keep track of what the previous bus was doing an 
hour ago, and yet one does not see clustering behavior on these low-frequency routes. The reason for this is that the 
instabilities predicted by Eq. [S] may take a long time to become noticeable, and normal bus routes tend to have a 
limited number of stops. Figure ITT1 shows the number of stops that a bus route has to cover before seeing a noticeable 
deviation in the initial homogeneous state. If the passenger rate is \i = 0.01 as suggested in the urban case above, 
then the route would have to have 130 stops to show a 1 minute deviation from the homogeneous state, and 225 stops 
to show a 5 minute deviation. If L = 0.5 km, these correspond to 65 km and 112 km, longer than your average bus 
route. The situation is even better when you consider that a suburban or rural bus route might have, not 1 person for 
every 5 minutes, but maybe 1 person every 25 minutes, so that [i = 0.002, and we can start having bus routes with 
600 or 1000 stops before the instability becomes noticeable. This is not to say that smaller routes remain perfectly 
on time, of course; just that the delays are unlikely to be due to the need to pick up extra passengers. Since buses 
will typically complete the route only to turn around and do it again, one might consider an entire day's run to be 
a single route, in which case instabilities may creep in late in the day. However, the introduction of a bus terminal, 
where buses wait until a specific time to leave for their next trip through the route, would have to be accounted for 
in this case. 

In this paper, we have considered the bus route model proposed in Ref. Q. We have added a simple way to deal 
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with negative time-headways (by replacing all negative Atj-'s with zeroes), and by doing so have been able to study 
the unstable modes of a homogeneous system of buses. We show that there are in fact three different phase diagrams 
(Fig El for this system, depending on our choice of parameters, and that in addition to the stable homogeneous 
state, there are three unstable modes which the system can fall into: the explosive mode, the slowed mode, and the 
underdamped mode. We then considered the application of this model to reahlife bus routes. We have shown that 
the passenger rate /i and the initial spacing between buses Ato are greatly restricted by practical considerations, and 
that to guarantee stability one needs to have bus drivers who are constantly tracking the bus in front of them, even 
when that bus is at its normal distance. Fortunately, this is only necessary for bus routes with very many stops; with 
fewer stops to make, a bus may be able to complete the route before any instabilities can become noticeable. 
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